Comparative genomics of macaques and integrated insights into genetic variation and population history.

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds.

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

A region of suppressed recombination misleads neoavian phylogenomics.

Genomes are typically mosaics of regions with different evolutionary histories. When speciation events are closely spaced in time, recombination makes the regions sharing the same history small, and the evolutionary history changes rapidly as we move along the genome. When examining rapid radiations such as the early diversification of Neoaves 66 Mya, typically no consistent history is observed across segments exceeding kilobases of the genome. Here, we report an exception. We found that a 21-Mb region in avian genomes, mapped to chicken chromosome 4, shows an extremely strong and discordance-free signal for a history different from that of the inferred species tree. Such a strong discordance-free signal, indicative of suppressed recombination across many millions of base pairs, is not observed elsewhere in the genome for any deep avian relationships. Although long regions with suppressed recombination have been documented in recently diverged species, our results pertain to relationships dating circa 65 Mya. We provide evidence that this strong signal may be due to an ancient rearrangement that blocked recombination and remained polymorphic for several million years prior to fixation. We show that the presence of this region has misled previous phylogenomic efforts with lower taxon sampling, showing the interplay between taxon and locus sampling. We predict that similar ancient rearrangements may confound phylogenetic analyses in other clades, pointing to a need for new analytical models that incorporate the possibility of such events.

Adaptive expansion of ERVK solo-LTRs is associated with Passeriformes speciation events.

Endogenous retroviruses (ERVs) are ancient retroviral remnants integrated in host genomes, and commonly deleted through unequal homologous recombination, leaving solitary long terminal repeats (solo-LTRs). This study, analysing the genomes of 362 bird species and their reptilian and mammalian outgroups, reveals an unusually higher level of solo-LTRs formation in birds, indicating evolutionary forces might have purged ERVs during evolution. Strikingly in the order Passeriformes, and especially the parvorder Passerida, endogenous retrovirus K (ERVK) solo-LTRs showed bursts of formation and recurrent accumulations coinciding with speciation events over past 22 million years. Moreover, our results indicate that the ongoing expansion of ERVK solo-LTRs in these bird species, marked by high transcriptional activity of ERVK retroviral genes in reproductive organs, caused variation of solo-LTRs between individual zebra finches. We experimentally demonstrated that cis-regulatory activity of recently evolved ERVK solo-LTRs may significantly increase the expression level of ITGA2 in the brain of zebra finches compared to chickens. These findings suggest that ERVK solo-LTRs expansion may introduce novel genomic sequences acting as cis-regulatory elements and contribute to adaptive evolution. Overall, our results underscore that the residual sequences of ancient retroviruses could influence the adaptive diversification of species by regulating host gene expression.

Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time.

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

Divergent Evolutionary Rates of Primate Brain Regions as Revealed by Genomics and Transcriptomics.

Although the primate brain contains numerous functionally distinct structures that have experienced diverse genetic changes during the course of evolution and development, these changes remain to be explored in detail. Here we utilize two classic metrics from evolutionary biology, the evolutionary rate index (ERI) and the transcriptome age index (TAI), to investigate the evolutionary alterations that have occurred in each area and developmental stage of the primate brain. We observed a higher evolutionary rate for those genes expressed in the non-cortical areas during primate evolution, particularly in human, with the highest rate of evolution being exhibited at brain developmental stages between late infancy and early childhood. Further, the transcriptome age of the non-cortical areas was lower than that of the cerebral cortex, with the youngest age apparent at brain developmental stages between late infancy and early childhood. Our exploration of the evolutionary patterns manifest in each brain area and developmental stage provides important reference points for further research into primate brain evolution.

Olympic gels formed through catenation of dsDNA rings regulated by topoisomerase II: A coarse-grained model.

Topological regulation of DNA by topoisomerases in cells is very crucial for life. We propose a coarse-grained model to study the catenation process of double-stranded DNA (dsDNA) rings regulated by topoisomerase II (TOP2) and provide a computational method to characterize the topological structures of the Olympic gels obtained. The function of TOP2 in the catenation of dsDNA rings is implicitly fulfilled by operating the length of a stretchable catch bond in the dsDNA ring. After the catenation reaction of initially noncatenated dsDNA rings in the solution, the Olympic gel is obtained and the interlocked topology of the dsDNA rings can be characterized by a computational method derived from the HOMFLY polynomial, based on which the catenation degree and the complexity of catenation are quantified. Detailed dependence of the catenation degree and the complexity of the catenated topology on key parameters, including the size of the transient broken gap and the duration time of the break on the dsDNA ring during operation by TOP2, the initial molar ratio of TOP2 to the dsDNA rings, and the reaction temperature, has been investigated.

Phylogenomic analyses unraveled the evolution of viral tolerance in bats.

Bats host a range of viruses, exhibiting a coevolution process with many virus genera and a special capacity for viral tolerance. Foley et al.1 performed phylogenomic analyses for 60 bat species, finding that swarming behavior might facilitate cross-species introgression and the spread of anti-virus immunity gene loci across species.

Color vision evolution in egg-laying mammals: insights from visual photoreceptors and daily activities of Australian echidnas.

Egg-laying mammals (monotremes) are considered "primitive" due to traits such as oviparity, cloaca, and incomplete homeothermy, all of which they share with reptiles. Two groups of monotremes, the terrestrial echidna (Tachyglossidae) and semiaquatic platypus (Ornithorhynchidae), have evolved highly divergent characters since their emergence in the Cenozoic era. These evolutionary differences, notably including distinct electrosensory and chemosensory systems, result from adaptations to species-specific habitat conditions. To date, very few studies have examined the visual adaptation of echidna and platypus. In the present study, we show that echidna and platypus have different light absorption spectra in their dichromatic visual sensory systems at the molecular level. We analyzed absorption spectra of monotreme color opsins, long-wavelength sensitive opsin (LWS) and short-wavelength sensitive opsin 2 (SWS2). The wavelength of maximum absorbance (λmax) in LWS was 570.2 in short-beaked echidna (Tachyglossus aculeatus) and 560.6 nm in platypus (Ornithorhynchus anatinus); in SWS2, λmax was 451.7 and 442.6 nm, respectively. Thus, the spectral range in echidna color vision is ~ 10 nm longer overall than in platypus. Natural selection analysis showed that the molecular evolution of monotreme color opsins is generally functionally conserved, suggesting that these taxa rely on species-specific color vision. In order to understand the usage of color vision in monotremes, we made 24-h behavioral observations of captive echidnas at warm temperatures and analyzed the resultant ethograms. Echidnas showed cathemeral activity and various behavioral repertoires such as feeding, traveling, digging, and self-grooming without light/dark environment selectivity. Halting (careful) behavior is more frequent in dark conditions, which suggests that echidnas may be more dependent on vision during the day and olfaction at night. Color vision functions have contributed to dynamic adaptations and dramatic ecological changes during the ~ 60 million years of divergent monotreme evolution. The ethogram of various day and night behaviors in captive echidnas also contributes information relevant to habitat conservation and animal welfare in this iconic species, which is locally endangered.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease.

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Engineering the activity and stability of ZIF-8(Zn/Co)@g-C3N4 nanocomposites and their synergistic action in converting atmospheric CO2 into cyclic carbonates.

The development of novel catalytic materials that integrate multifunctional sites has significant implications for expanding the utilization of CO2 resources. However, simultaneously achieving high activity and stability remains a formidable challenge. In this study, a series of ZIF-8(Zn/Co)@g-C3N4 nanocomposites were prepared by employing a thermo-physical compounding strategy that involved the combination of nitrogen-rich graphitic carbon nitride (g-C3N4) nanosheets with ZIF-8(ZnCo). The influences of different compositions of g-C3N4 and ZIF-8(Zn/Co) on the catalyst structure were systematically investigated. Subsequently, the catalytic activities of these nanocomposites towards the cycloaddition reaction between CO2 and epoxide were examined under different conditions. The presence of abundant Lewis base sites in g-C3N4 facilitates CO2 activation, while multiple Lewis acid sites in ZIF-8(Zn/Co) enable efficient epoxide activation. By working synergistically with a co-catalyst, tetrabutylammonium bromide (TBAB), CO2 and epoxides can be efficiently reacted to synthesize the corresponding cyclic carbonates under mild or even atmospheric pressure conditions. The catalytic reaction conditions were optimized, and both the catalyst's recycling performance and the scope of epoxides with various substituents were investigated. The integration of g-C3N4 and ZIF-8(Zn/Co) endows the catalytic material with exceptional structural stability and remarkable catalytic activity, thereby providing a new platform for highly efficient CO2 conversion.

High-Performance Ethylene Glycol Sensor Based on Imine Covalent Organic Frameworks.

The colorless and odorless ethylene glycol is prone to unknowingly causing poisoning, making preventive monitoring of ethylene glycol necessary. In this paper, scandium (III) trifluoromethanesulfonate was used as a catalyst to successfully prepare covalent organic framework (COF) nanospheres linked by imines at room temperature. The COF nanospheres were characterized by XRD, SEM, TEM, FT-IR, UV-Vis and BET. The results show that COF nanospheres have rough surfaces and a large number of mesoporous structures, which greatly increase the active sites on the surface of the sensing material and enhance the gas sensing performance. The sensing results showed that the prepared imine-conjugated COF nanospheres exhibited a good response-recovery ability for 10 consecutive response-recovery cycles for ethylene glycol at room temperature and had a theoretical detection limit of 40 ppb. In addition, the responses of COF nanospheres to nearly 20 interfering gases, including HCl, HNO3, phenol, formaldehyde and aniline, are relatively low compared to the response to ethylene glycol, indicating that the COF nanospheres have high selectivity towards ethylene glycol. The COF nanospheres show good sensitivity and selectivity for the detection of ethylene glycol, which should be attributed to the large specific surface area, hydrogen bonding interactions, and high defects. This work provides an effective method for the detection of ethylene glycol and expands the application field of COF materials.

Genetic Load and Adaptive Potential of a Recovered Avian Species that Narrowly Avoided Extinction.

High genetic diversity is a good predictor of long-term population viability, yet some species persevere despite having low genetic diversity. Here we study the genomic erosion of the Seychelles paradise flycatcher (Terpsiphone corvina), a species that narrowly avoided extinction after having declined to 28 individuals in the 1960s. The species recovered unassisted to over 250 individuals in the 1990s and was downlisted from Critically Endangered to Vulnerable in the International Union for the Conservation of Nature Red List in 2020. By comparing historical, prebottleneck (130+ years old) and modern genomes, we uncovered a 10-fold loss of genetic diversity. Highly deleterious mutations were partly purged during the bottleneck, but mildly deleterious mutations accumulated. The genome shows signs of historical inbreeding during the bottleneck in the 1960s, but low levels of recent inbreeding after demographic recovery. Computer simulations suggest that the species long-term small Ne reduced the masked genetic load and made the species more resilient to inbreeding and extinction. However, the reduction in genetic diversity due to the chronically small Ne and the severe bottleneck is likely to have reduced the species adaptive potential to face environmental change, which together with a higher load, compromises its long-term population viability. Thus, small ancestral Ne offers short-term bottleneck resilience but hampers long-term adaptability to environmental shifts. In light of rapid global rates of population decline, our work shows that species can continue to suffer the effect of their decline even after recovery, highlighting the importance of considering genomic erosion and computer modeling in conservation assessments.

Dissecting aneuploidy phenotypes by constructing Sc2.0 chromosome VII and SCRaMbLEing synthetic disomic yeast.

Aneuploidy compromises genomic stability, often leading to embryo inviability, and is frequently associated with tumorigenesis and aging. Different aneuploid chromosome stoichiometries lead to distinct transcriptomic and phenotypic changes, making it helpful to study aneuploidy in tightly controlled genetic backgrounds. By deploying the engineered SCRaMbLE (synthetic chromosome rearrangement and modification by loxP-mediated evolution) system to the newly synthesized megabase Sc2.0 chromosome VII (synVII), we constructed a synthetic disomic yeast and screened hundreds of SCRaMbLEd derivatives with diverse chromosomal rearrangements. Phenotypic characterization and multi-omics analysis revealed that fitness defects associated with aneuploidy could be restored by (1) removing most of the chromosome content or (2) modifying specific regions in the duplicated chromosome. These findings indicate that both chromosome copy number and specific chromosomal regions contribute to the aneuploidy-related phenotypes, and the synthetic chromosome resource opens new paradigms in studying aneuploidy.

Validation of key sponge symbiont pathways using genome-centric metatranscriptomics.

The sponge microbiome underpins host function through provision and recycling of essential nutrients in a nutrient poor environment. Genomic data suggest that carbohydrate degradation, carbon fixation, nitrogen metabolism, sulphur metabolism and supplementation of B-vitamins are central microbial functions. However, validation beyond the genomic potential of sponge symbiont pathways is rarely explored. To evaluate metagenomic predictions, we sequenced the metagenomes and metatranscriptomes of three common coral reef sponges: Ircinia ramosa, Ircinia microconulosa and Phyllospongia foliascens. Multiple carbohydrate active enzymes were expressed by Poribacteria, Bacteroidota and Cyanobacteria symbionts, suggesting these lineages have a central role in assimilating dissolved organic matter. Expression of entire pathways for carbon fixation and multiple sulphur compound transformations were observed in all sponges. Gene expression for anaerobic nitrogen metabolism (denitrification and nitrate reduction) were more common than aerobic metabolism (nitrification), where only the I. ramosa microbiome expressed the nitrification pathway. Finally, while expression of the biosynthetic pathways for B-vitamins was common, the expression of additional transporter genes was far more limited. Overall, we highlight consistencies and disparities between metagenomic and metatranscriptomic results when inferring microbial activity, while uncovering new microbial taxa that contribute to the health of their sponge host via nutrient exchange.

Historic Sampling of a Vanishing Beast: Population Structure and Diversity in the Black Rhinoceros.

The black rhinoceros (Diceros bicornis L.) is a critically endangered species historically distributed across sub-Saharan Africa. Hunting and habitat disturbance have diminished both its numbers and distribution since the 19th century, but a poaching crisis in the late 20th century drove them to the brink of extinction. Genetic and genomic assessments can greatly increase our knowledge of the species and inform management strategies. However, when a species has been severely reduced, with the extirpation and artificial admixture of several populations, it is extremely challenging to obtain an accurate understanding of historic population structure and evolutionary history from extant samples. Therefore, we generated and analyzed whole genomes from 63 black rhinoceros museum specimens collected between 1775 and 1981. Results showed that the black rhinoceros could be genetically structured into six major historic populations (Central Africa, East Africa, Northwestern Africa, Northeastern Africa, Ruvuma, and Southern Africa) within which were nested four further subpopulations (Maasailand, southwestern, eastern rift, and northern rift), largely mirroring geography, with a punctuated north-south cline. However, we detected varying degrees of admixture among groups and found that several geographical barriers, most prominently the Zambezi River, drove population discontinuities. Genomic diversity was high in the middle of the range and decayed toward the periphery. This comprehensive historic portrait also allowed us to ascertain the ancestry of 20 resequenced genomes from extant populations. Lastly, using insights gained from this unique temporal data set, we suggest management strategies, some of which require urgent implementation, for the conservation of the remaining black rhinoceros diversity.

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Tuning Lignite Structure via Hydromodification To Promote the Formation of Coal-Based CNTs: Exploration for the Carbon Source of CNTs.

Although the preparation of coal-based carbon nanotubes (CNTs) has been realized in many studies, the relationship between carbon source structure of coal and CNT growth has not been studied in depth. In this study, we used lignite and KOH as raw material and catalyst and tuned lignite structure via hydrothermal modification to promote the formation of CNTs during catalytic pyrolysis. The main carbon source of CNTs was explored from the change of coal structure and pyrolysis characteristics. The results indicate that the CNT yield of lignite pyrolysis products is only 2.39%, but the CNT yield increases significantly after lignite was hydrothermally modified in a subcritical water-CO system. The graphitization degree, the order degree, and CNT content increase continuously with the increase in modification temperature, and C-M340 has the highest CNT content of 9.41%. Hydromodification promotes the rearrangement of aromatic carbon structures to generate more condensed aromatic rings linked by short aliphatic chains and aromatic ether bonds. The variation of these structures correlates well with the formation of CNTs and leads to the change in the carbon source components released during coal pyrolysis. Compared to lignite, modified coal releases more aromatic compounds, especially polycyclic aromatic hydrocarbons with ≥3 rings and phenols during catalytic pyrolysis, which is conducive to the transformation into carbon clusters and provides carbon sources for CNT growth. In addition, modified coal releases a slightly more carbon-containing gas (CH4 and CO) than lignite, which has a limited effect on the growth of CNTs. This study provides a novel and efficient method for enhancing the growth of CNTs by a molecular tailoring strategy of coal.

Integrative Omics Reveals Rapidly Evolving Regulatory Sequences Driving Primate Brain Evolution.

Although the continual expansion of the brain during primate evolution accounts for our enhanced cognitive capabilities, the drivers of brain evolution have scarcely been explored in these ancestral nodes. Here, we performed large-scale comparative genomic, transcriptomic, and epigenomic analyses to investigate the evolutionary alterations acquired by brain genes and provide comprehensive listings of innovatory genetic elements along the evolutionary path from ancestral primates to human. The regulatory sequences associated with brain-expressed genes experienced rapid change, particularly in the ancestor of the Simiiformes. Extensive comparisons of single-cell and bulk transcriptomic data between primate and nonprimate brains revealed that these regulatory sequences may drive the high expression of certain genes in primate brains. Employing in utero electroporation into mouse embryonic cortex, we show that the primate-specific brain-biased gene BMP7 was recruited, probably in the ancestor of the Simiiformes, to regulate neuronal proliferation in the primate ventricular zone. Our study provides a comprehensive listing of genes and regulatory changes along the brain evolution lineage of ancestral primates leading to human. These data should be invaluable for future functional studies that will deepen our understanding not only of the genetic basis of human brain evolution but also of inherited disease.

Eighty million years of rapid evolution of the primate Y chromosome.

The Y chromosome usually plays a critical role in determining male sex and comprises sequence classes that have experienced unique evolutionary trajectories. Here we generated 19 new primate sex chromosome assemblies, analysed them with 10 existing assemblies and report rapid evolution of the Y chromosome across primates. The pseudoautosomal boundary has shifted at least six times during primate evolution, leading to the formation of a Simiiformes-specific evolutionary stratum and to the independent start of young strata in Catarrhini and Platyrrhini. Different primate lineages experienced different rates of gene loss and structural and chromatin change on their Y chromosomes. Selection on several Y-linked genes has contributed to the evolution of male developmental traits across the primates. Additionally, lineage-specific expansions of ampliconic regions have further increased the diversification of the structure and gene composition of the Y chromosome. Overall, our comprehensive analysis has broadened our knowledge of the evolution of the primate Y chromosome.